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It is well known that at the gel-liquid phase transition temperature a lipid bilayer membrane 
exhibits an increased ion permeability. We analyze the quantized currents in which the increased 
permeability presents itself. The open time histogram shows a "-3/2" power law which implies an 
open-closed transition rate that decreases like k{t) oc t~^ as time evolves. We propose a "pore 
freezing" model to explain the observations. We discuss how this model also leads to the 1//" noise 
that is commonly observed in currents across biological and artificial membranes. 

PACS numbers: 05.70.Fh, 87.10.Mn, 87.16.dp 



I. INTRODUCTION 

Lipid bilayer membranes undergo melting transitions 
at temperatures that are generally close to physiological 
temperatures p. In the fluid phase there is free lateral 
diffusion of the lipids with a diffusion constant in the 
lim^/s regime. In the solid or gel phase the lipids are 
arranged in a more rigid two-dimensional lattice. Near 
the phase transition temperature solid lipid domains drift 
in the liquid embedding. 

The solid-to-liquid transition involves significant 
changes in the volume and in the surface area. This leads 
to increased volume compressibility and area compress- 
ibility at the phase transition [Tj. Consequently, Brow- 
nian noise causes relatively large fluctuations in volume 
and area at the point of phase transition. Monte-Carlo 
simulations have affirmed that fluctuations are indeed 
particularly large at the solid-liquid interface near solid 
domains in the liquid membrane [5j . 

Already in 1973 it was reported that the sodium per- 
meability of an artificial lipid bilayer membrane peaked 
sharply at the melting temperature P]- Given the large 
fluctuations at the solid-liquid interface this result is no 
surprise. What was truly surprising was the discovery 
in 1980 by Antonov et al that this increased permeabil- 
ity comes in the form of quantized currents [Ij. For 
a relatively small transmembrane voltage (~100 mV) 
"channels" of a fixed conductance in the picosiemens to 
nanosiemens regime appeared to open and close. The 
behavior strongly reminds of that of ion channel proteins 
[5]. It is startling that just the lipid domains by them- 
selves can actually exhibit the same behavior as is exhib- 
ited by the specialized complex proteins that ion channels 
are. The results of Antonov et al have been reproduced 
many times and for many different kinds of membranes 
(see references in [T]). However, the quantized currents 
that appear to turn on and off have, so far, defied ex- 
planation. In this article we study the statistics of the 
open and closed times. On the basis of the analysis of 
experimental data we propose an explanation. 

Electroporation, i.e. the "punching" of holes in a lipid 
bilayer membrane with an electric field, has been studied 
extensively by theoreticians and experimentalists [B]. In 



our experiments we followed a procedure that was first 
described by Miiller et al [7|: a black lipid membrane 
(BLM) is assembled in a small aperture (80 /im diame- 
ter) in a thin (25 ^m) electrically insulating teflon film. 
The compartments on both sides of the film contain a 
150 mM KCl solution. Across the film there is a constant 
electric potential and the resulting transmembrane cur- 
rent is measured. For further details the reader is referred 
to Refs. '8 and [IT. Our bilayer membrane was made up 
of dioleoyl-phosphatidylcholine (DOPC) and dipalmitoyl 
phosphatidylcholine (DPPC) in a 2:1 ratio. We picked 
these phospholipids and this ratio because they give rise 
to a phase transition at room temperature, i.e. 19°C. For 
reasons explained later in this article, the membrane was 
also made to contain 15.9 mol% octanol. Typical results 
are depicted in Fig. 1. 
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FIG. 1. A typical result of the BLM experiment described in 
the text (see also Refs. [5] and [3]). This current trace was 
obtained at a 210 mV transmembrane voltage. Measurements 
were done at a few degrees above the melting temperature. 
We sampled data for 30 seconds. A 300 Hz filter was applied. 
The resulting open times are collected in a histogram (Fig. 
2a). 



II. DATA AND ANALYSIS 

Figure 1 does indeed look very much like an ion channel 
that is fluctuating between one open and one closed state. 
We put the cut-off at at /™ — 10 pA and associated 
a transmembrane current of Im < Im with the closed 
state and a transmembrane current of Im > Im with the 
open state. The results that we found on the statistics 
of open and closed time intervals appeared robust under 
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FIG. 2. (a) A histogram on a log- log scale showing the rate 
of occurrence of the different pore lifetimes. In order to cover 
the wide range of lifetimes, different bin sizes where combined 
in one graph following a procedure described in Ref . nTT . For 
lifetimes larger than about 3 x 10~^ s the data arc well fitted 
with the solid line that represents the theoretically inferred 
f{t) = t~'^'^ /{2yfK\). (b) The derived open-to-closed tran- 
sition rate for an electropore. These data appear to be well 
fitted with the solid line k{t) = l/{2t), which is the power law 
that the theory predicts. 

small variations (about ±1 pA) of /™. We took data for 
30 seconds. We checked that there was no drift of the 
open-closed statistics in the course of the analyzed time 
interval. 

When putting together an open time histogram like 
Fig. 2a, the bin size choice generally constitutes a com- 
promise. For a small bin size, the number of events in 
each bin may become very small and, especially for large 
open times, the statistical variations from bin to bin will 
overwhelm any general trend. For larger bin size, one 
can lose a lot of information in the region of the shorter 
open times. In Ref. [TD] a method is described to join 
data for different bin sizes into one histogram. In Ref. 
[llj an actual computer program is provided to construct 
such a histogram. That program was configured to work 
with our data. At each bin width, twenty bins were al- 
located to count occurrences. The smallest bin size, i^, 
was twice the minimum open time and at each cycle, the 
bin size increased as ^i+i — 1.1 xti. In Ref. [TT] a factor 
2 was used in place of 1.1. We took the smaller number 
so as to generate more data points over a smaller range 
of time. In order to fit the histograms with different bin 
sizes together, the counts in each bin were normalized 
by 1 1 {tiNtotai) , where Ntotai is the total number of open 
times. The first bin must be excluded because it will 
contain all of the the times unresolved at time ti [11] . As 
soon as a bin is encountered that contains no counts the 
loop breaks to the next ti. Eventually, the bins become 
large enough that the counts become too sparse or all of 
the counts are in the first bin and no more information 
can be obtained. Figure 2a shows how, over more than 
two decades, the open time histogram is well fitted with 
a power law. 

When normalized, the histogram is a probability den- 
sity function of open times f{t). If we let P{t) denote the 
probability that the pore stays open till time t or longer. 



then f{t) = —P{t). The distribution f{t) is related to 
the transition rate k{t) out of the open state in the fol- 
lowing way: k{t) = -f{t)/P{t) = -J^ In P(t). Reference 
[TTj also shows how to calculate the transition rate k{t) 
from the open time histogram data. Kinetic rate con- 
stants were calculated at each bin size by fitting the 2"**, 
3'''^, and 4*'' bins to a single exponential if those bins 
contained any counts. The result is depicted in Fig. 2b. 

The pure lipid bilayer has a heat capacity maximum 
at 19° C. The phase transition occurs in a range around 
this temperature. At 19° C we regularly recorded what 
looked like two channels that were open at the same time. 
The analysis of such observations is problematic: when 
one of two simultaneously open channels closes, there is 
no way to know which one of the two it is. The result- 
ing ambiguity makes an accurate open time histogram 
impossible. We also did a number of experiments with 
different concentrations of octanol dissolved in the mem- 
brane. The anesthetic octanol is a very nonpolar sub- 
stance. Much of this substance will dissolve among the 
lipid tails of the phospholipids. This will bring down 
the melting temperature, , of the membrane through 
the same mechanism by which the addition of salt brings 
down the freezing temperature of water [T^ . By dissolv- 
ing the appropriate amount of octanol in the membrane, 
we can achieve any desired shift of . Increasing the dif- 
ference between the ambient temperature and makes 
channel openings correspondingly more rare. The results 
that are depicted in our Figs. 1 and 2 were obtained with 
an octanol concentration in the membrane of 15.9 mol%. 
In Fig. 3 of Ref. [5] it can be seen that 15.9 mol% octanol 
leads to a melting temperature that is about 5°C below 
the temperature at which experiments were performed. 
With 15.9 mol% octanol double openings constituted a 
negligible fraction of the total amount of openings. Sim- 
ilar results were obtained for 7.9 mol% octanol [3 [9]. 

Ion channel proteins with a closed state (C) and an 
open state (O) have commonly been modeled as two- 
state molecules with Markov transitions connecting the 
two states. The simple C^O scheme with transition 
rates that do not depend on time leads to exponentially 
distributed lifetimes in both the open and closed state. 
However, nonexponential distributions like our apparent 
power-law (Fig. 2a) have been often encountered. Elec- 
trophysiologists have commonly explained nonexponen- 
tial distributions of closed and open times with kinetic 
schemes that contain more than two states |13j . 

An alternative approach has been to model a transition 
with a time dependent transition rate, e.g. k{t) cx 
[lOj . There is ample biomolecular justification for such 
a decreasing transition rate when we are dealing with 
the complicated ion channel proteins. These are protein 
complexes with a behavior that is often not adequately 
described with states and rates or with a straightforward 
ID reaction coordinate. Even a relatively simple protein 
like myoglobin turns out to be better described in terms 
of diffusive motion in a many-dimensional conformational 
space [14 . With the power law k{t) oc we imagine 
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a channel that, after crossing the activation barrier to 
the other state, diffuses away from that activation bar- 
rier. The transition rate back to the original state then 
decreases and the channel thus "stabilizes in its state" as 
time progresses. 

In 1988 Millhauser, Salpetcr, and Oswald proposed the 
following kinetic scheme to explain the Fig. 2 type statis- 
tics that had, at that time, been found for closed times 
of several ion channel proteins |15j : 
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This model applies to our case in the following way. State 
"0" denotes the initial absence of the electropore. The 
pore forms amid high fluctuations in the liquid area near 
a solid-liquid interface. The states 1 to denote open 
states of the pore. We propose that going from 1 to 
corresponds to the freezing of an increasing amount of 
lipids in and around the pore's lining. Upon transiting 
to state 2, a pore can only close again after it goes back 
to state 1. The states 1 to TV all correspond to the same 
pore size and conductance. But as the state of the pore 
moves from 1 to A'^, the pore becomes embedded in an 
ever larger solid region. For the pore to close again, it 
is necessary that the "frozen" surroundings melt again. 
Such melting can only occur from the boundaries inward 
towards the pore. In terms of the kinetic scheme Q , the 
decreasing k{t) for the closing transition comes about as 
a significant fraction of the probability diffuses away to- 
wards A^. Though the individual transitions in the above 
kinetic scheme are Markovian, there is a time-dependent 
rate k{t) from the 1-to-A^ set of states back to the closed 
0-state. 

More than a decade after the above kinetic scheme ([T]) 
was first proposed, Goychuk and Hanggi published a se- 
ries of papers where the discrete diffusion in terms of 
states and rates was replaced by a continuous time ran- 
dom walk [16] . They found that the closed time statistics 
were the same in this case as what Millhauser, Salpeter, 
and Oswald had derived. Even for an energy landscape 
that is not flat, but has some minor undulations they 
found the same power laws to hold. The model that 
was worked out by Millhauser, Salpeter, and Oswald and 
by Goychuk and Hanggi (MSOGH) was intended to de- 
scribe diffusion in a protein's conformational space. Such 
a space is an abstraction. But in our case the MSOGH 
model actually describes the much less abstract growth 
of a solid domain. Simulations presented in [3] show that 
even at a few degrees above T„i, solid domains form in 
the membrane. 

We assume that a pore will most readily form at a 
solid-liquid interface where area and volume fluctuations 
are high. A pore can then next stabilize if the involved 
phospholipids freeze. It will stabilize further if phospho- 
lipids around the pore also freeze. It is this freezing that 
can make a pore long lived. This model assumption is 
consistent with the observation that pores in a purely liq- 
uid membrane have much higher activation barriers for 
their formation (more than 250 mV is generally required 



for observable pore formation) and close within millisec- 
onds after the transmembrane voltage is brought back to 
physiological levels [T71 [T5| . 

It turns out that the kinetic scheme ([T]) can be analyt- 
ically solved for the case in which all rate constants are 
identical [H]. This is exactly the case we have at hand 
between states 1 and A^, as our diffusion is isotropic and 
has no directional bias. Going from state 1 into the di- 
rection of state A^ is just the attachment of more lipids 
to the solid raft. For noise suppression we used a 300 Hz 
fllter as we analyzed the data. This means that, in effect, 
we are taking a snapshot of the system about every three 
milliseconds. It is therefore that Fig. 2 shows a break- 
down of the power law at about 3 ms. Because of the 300 
Hz cut-off there is no way of knowing how far the power 
law extends for the actual physical process. We therefore 
have to take A w 300 s~^ as the rate for the "0 1"- 
transition in the kinetic scheme ([T]). We can model the 
fluctuating size of the solid domain as a diffusion in con- 
formational space and, with the aforementioned results 
of Goychuk and Hanggi [16], we take all the rates in the 
kinetic scheme to be the same A ss 300 s~^. 

When all the rates are normalized, the master equation 
that is associated with the kinetic scheme is: 

-2pi +P2 

Pn{t)^ Pn-i - 2p„ -I- Pn+1 for 1 < n < A^ (2) 

PN{t)^ PN-l -~PN ■ 

We have for the aforementioned P{t): 

Pit) = S,lip„(i) . (3) 

It is obvious from the kinetic scheme ([T]) that f[t) — 
—dP{t)/dt = pi{t). We start out with pi{0) — 1 and 
p„(0) = for n > 2. This imphes P{0) = 1. For A^ -> oo 
the analytic solution of Eq. (2) is: 



Pn{t) = exp[-2t] (In-im - In+l{2t)) , 



(4) 



where /„(i) represents the modified Bessel function of 
order n. An expansion for t » ^n? leads to: 



Pr, 



(5) 



For the open-closed transitions, all that matters is pi{t). 
So for < > 1 we have a good approximation with: 
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All of these results are also shown in [T5] . 

We have to unscale the normalized transition rates in 
the kinetic scheme before we can relate the data in Fig. 
2 to the theory. Redimensionalization turns the t 1 
condition that goes with Eq. (|6| into t > 3 ms. Unsealing 



Eq. (|6|) leads to f{t) w r-^'^ /{2\/'kX). As k{t) x t is 
dimensionless the k{t) « l/(2t) remains unaffected by 
the unsealing. Figure 2 shows how the theory provides 
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FIG. 3. A pore in a lipid bilayer involves a rearrangement 
of a number phospholipids. The figure is adapted from R.A. 
Bockmann, B.L. de Groot, S. Kakorin, E. Neumann, and H. 
Grubmiiller, Biophys. J. 95, 1837 (2008). The polar head- 
groups of the phospholipids make up the lining of the pore and 
make the pore permeable to water and small ions. Shown is 
a schematic idea of (a) the narrowest part of the pore viewed 
from the direction perpendicular to the membrane and (b) a 
cross sectional view in the plane of the membrane through 
the center of the pore. 



an almost perfect fit to the data. The kinetic scheme ([T]) 
and the ensuing theory correctly predict both the slope 
and the height of the data points. 

The current traces (Fig. 1) make clear that the conduc- 
tance of a single pore has a preferred level. Histograms 
presented in Refs [5] and [3] make this assertion more rig- 
orous. Figure 3 gives a schematic idea of the architecture 
of a pore in a lipid bilayer. If we assume the pore to be 
cylindrical and electrolyte filled, then the conductance of 
70 pS that we have at 15.9 mol% octanol corresponds to 
a pore of about 0.35 nm radius. Such a radius is indeed 
similar to that of many protein channels. A single head- 
group covers a surface area of about 0.6 nm^ in the fluid 
state and of about 0.5 nm^ in the gel state. Assuming 
pores again to be roughly cylindrical and the pore length 
to be 5 nm, this means that a pore of a 0.35 nm radius 
involves about 20 phospholipids. 

Considerations involving standard electropore theory 
make clear why this quantization can occur. Forming a 
pore involves a rearrangement of phospholipids and an 
activation barrier has to be crossed to bring about such 
rearrangement |19j . Obviously, a minimum pore radius 
ro of a few tenths of nanometers is required for a pore 
to be permeable to water and ions. For phosphate head- 
groups to keep facing the water even in the pore's inte- 
rior, the phospholipids have to form a very curved edge 
(see Fig. 3). The energy required to create the edge of 
a pore of radius r is 2Trjr, where 7 is the so-called line 
tension [3D]. For small pores we can approximate the 
energy with just the linear terms in the radius r, i.e., 
E{r) = (27r7 — eoei^V'^)r. Here V is the transmembrane 
voltage and e^, is the relative dielectric permittivity of 
water. The term —eQEwV^ describes the Maxwell stress 
due to the inhomogeneity of the electric field that the 
conducting aqueous pore causes [3T]. This stress is to- 
wards further opening of the pore. Electroporation has 
commonly been studied with voltages over 300 mV. Such 
voltages, after all, are required to permcabilize a lipid bi- 
layer in its liquid state. But for the smaller voltages that 
we work with the energy E increases very rapidly with r. 



So, in our case, thermal fluctuations (with a free energy 
of the order of ksT) will not be sufficiently strong to 
drive the pore to radii that are significantly larger than 
tq. Once frozen a pore will, of course, not change its 
radius. 

We also made histograms for the closed times. If every 
apparent pore opening in Fig. 1 were to stem from the 
formation of a new pore, then we would have a constant 
pore formation rate. Such a constant closed-to-open rate 
would lead to an exponential distribution of closed times. 
However, also for the closed times we observe a power law 
and an /(t) and k(t) that appear to follow Eq. This 
non-Markov behavior indicates that a poreless membrane 
carries a memory for how long ago the most recent pore 
closed. This puzzling phenomenon could be explained 
with the realization that the lining of a pore consists of 
about 20 phospholipids and that it is possible for the 
pore lining to partially melt. In that case some individ- 
ual lipids in the lining could possibly get unstuck, move 
more freely, and clog up the pore. Meanwhile the main 
pore architecture remains in place and the pore can open 
again and refreeze. When several molecules are able to 
move independently into and out of the pore's interior 
we have a mechanism that can be described by the ki- 
netic scheme ([T]) and an MSOGH model. The apparent 
"-3/2" -power law for the closed times may be due to a 
significant fraction of the observed openings actually be- 
ing such "re-openings." Molecules or clusters of molecules 
from the lining of the pore that are moving in and out of 
the pore's interior have also been offered as an explana- 
tion for the 1//" noise that has been observed in other 
biological and artificial channels in membranes 

Figure 4 shows how the power spectrum of our record- 
ing (Fig. 1) follows a power law in the about two decades 
between the millisecond and the second regime where 
power laws applied in Fig. 2. It can be derived that the 
power spectral density follows 1//^ when the residence 
times in both the open and closed state follow a distri- 
bution f{t) cx ^-(1+^) where < £> < 1 [23]. Our power 
spectrum indeed displays the predicted -1/2 slope to a 
good approximation. Noise with a power spectrum that 
drops off like 1//", where 0.5 < a < 1.5, is commonly 
characterized as "1// noise." [l^. Already in 1966 it was 
found that the membrane voltage at the node of Ranvier 
of a live nerve cell displays 1//" noise (with an a very 
close to unity) in the 10 to 1000 Hz regime [IS]. Re- 
cently, these data have been more accurately rerecorded 
and it has been claimed that the apparent power law 
comes about as a sum of Lorentzian contributions of in- 
dividual types of ion channel proteins [55] . However, the 
phenomenon discussed in this article may be partly be- 
hind this observation. 



III. DISCUSSION 

A 1/f noise spectrum has also been found when an 
electropore is current clamped [27]. But it should be re- 
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FIG. 4. The power spectral density as a function of the 
frequency for the signal depicted in Fig. 1. A straight line 
of a slope -1/2 is added for reference. Power law behavior, 
S{f) oc f~^^^, is apparent in the significant frequency range. 
It is because of the logarithmic scale that data points get more 
dense towards the right end of the graph. 



alized that the underlying physics is completely different 
for a current clamped electropore. When the current, /, 
through an electropore in a membrane is kept constant, a 
negative feedback mechanism is at work that keeps that 
pore open. This can be imagined as follows. Suppose 
that, through a Brownian fluctuation, a pore narrows. 



This would increase the pore's electrical resistance R. 
With V = IR, this would increase the transmembrane 
voltage V. Through the aforementioned Maxwell stress, 
such larger transmembrane voltage increases the force 
that widens the pore [5T]. It is possible to keep pores 
open for hours with a current clamp setup \27\ I28j . How 
this negative feedback mechanism precisely leads to 1/f 
noise is unclear, but 1/f noise is commonly observed in 
current clamp setups. Ordinary solid resistors have been 
shown to exhibit 1/f noise over six decades when current 
clamped [H]. It should be reahzed that our 1/f noise 
does not originate in such fluctuations around the open 
level. Our experiments are done under voltage clamp and 
our noise spectrum arises from the openings and closings 
of the pore, such as the ones that are are shown in Fig. 
1. 

We have presented a description and an explanation 
of ion channel like behavior of electropores in a lipid bi- 
layer membrane. The possibility has to be seriously con- 
sidered that some of the behavior that is traditionally 
attributed to ion channels may actually originate in the 
lipid membrane. Many important ion channels exhibit 
open times in the millisecond regime. When we extrap- 
olate the power laws found in Figs. 2 and 4 to the mil- 
lisecond regime, we see that frequent openings can occur 
there. 
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